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METHOD FOR CONTROLLING COMBUSTION IN AN INTERNAL 
COMBUSTION ENGINE AND PREDICTING PERFORMANCE AND 

EMISSIONS 

Field of the Invention 

The present invention relates to a method for controlling 
combustion, performance and emissions of a direct injection internal 
combustion engine and predicting the behavior of a direct injection 
5 internal combustion engine. An embodiment of the subject method is 

particularly useful for gaseous-fuelled direct injection internal combustion 
engines. 

Background of the Invention 

1 0 Accurate prediction of the behavior of intemal combustion engines 

without engine testing can be valuable. First, predicting how an engine 
will perform imder a set of conditions can reveal valuable infomiation 
without having to incur the cost and time involved in testing an actual 
engine. Second, expensive prototype refinement is not unusual when 

1 5 testing a new engine design. Predicting behavior of a concept without 
having to manage prototype development can be a helpful step in engine 
research accommodated by methods of predicting engine behavior. 
Accurate engine behavior prediction helps overcome these problems by 
initially testing new ideas for a new engine design or simulating an 

2 0 existing engine architecture prior to any investment in costly prototype 
parts or time consuming test cell work. Finally, methods of predicting 
engine behavior can allow for accelerated testing of a potential system 
without limitation of manpower and facility availability. While not a 
substitute for actual testing, accurate prediction can eliminate many 

25 otherwise reasonable ideas before an expensive test engine is built or 
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before valuable test cell time is used to test an existing engine. 

Additionally, as computer processing power continues to improve, 
open loop controls using engine prediction data can eliminate both closed 
loop control methods and reliance on expensive monitoring equipment and 
5 the associated costly electrical systems, both of which are prone to wear. 

Finally, as emissions controls become increasingly important for 
environmental reasons, it is valuable to provide a controller for closed 
loop on->board control adjusting emissions in response to predicted 
behavior for a set of user demanded operational conditions. Also, as is the 
1 0 case in engine design noted above, a predictive system can help to filter 
out new designs that fail to meet acceptable emission levels. 

Prior art in this area has provided, in general, systems that predict 
aspects of engine operation but fail to satisfactorily integrate these to 
provide a predictive engine. By way of example, combustion models have 
1 5 been developed that predict combustion within a combustion chamber on 
the basis of operator inputted values. The same is true of emissions 
models, injector models and engine cycle models. 

Examples of engine combustion models development are 
summarized in Heywood, J.B., Internal Combustion Engine Fundamentals 

2 0 (New York: McGraw-Hill, 1988). The combustion model is typically 

integrated into the engine cycle model. Often accompanying these 
combustion models are emission formation models that provide an 
indication of NOx and particulate matter emission. Engine cycle 
simulation is used to compute thermodynamics and gas dynamics in the 
25 engine. GT-Power^'^ is an example of an engine cycle process modeling 
module. For engine cycle prediction, several methods have been provided 
in such systems as StarCD™, Fluent''^'^, FIRE'^^ and VECTIS^*^. 
Hydsim^'" software is an example of a prior art injector simulation model. 
The software can compute injection delay, injection rate shape by 

3 0 modeling details of the injector, including its hydraulic system, 

mechanical system and electronic system. 



Each of these components has been developed, to a greater or 
lesser extent, but, in general, remain independent models of one aspect of 
engine behavior. By way of example, consider A. Chow et al., 
Thermodynamic modeling of complete engine systems - a review (London: 
Professional Engineering Publishing, 1999) and T. Morel et al., ''Virtual 
Engine/Powertrain/Vehicle" Simulation Tool Solves Complex Interacting 
System Issues (Warrendale, PA: Society of Automotive Engineers 
Intemational, 2003) 2003-01-0372. In each case, examples of models used 
to describe aspects of an engine are provided. 

Chow discusses the system integration of models which determine 
engine behavior for a specified rate of combustion. These models 
quantitatively represent heat transfer, fluid flow, turbo-machinery, and 
emissions production. As one exanq>le of a system integration, Morel 
describes software which, for a specified rate of combustion, provides a 
system integration of separate models for heat transfer, fluid flow, turbo- 
machinery, and hydraulic, electrical, thermal, mechanical, and controls 
elements of the engine system. Other prediction methods discussed by 
Morel include components with turbocharger, supercharger and 
intercooler support; co-simulation options; transient pipe temperature 
analysis; variable pipe length; multi-user database and dynamic valve 
motion. The methods taught emphasize the relationship of the engine 
valves with other components. 

Generally, the integrations discussed do not explicitly mclude the 
interaction between fuel injection behavior and the rate of combustion. 
Full engine simulation integration benefits if the method allows for the 
interdependence of injector performance, combustion pressure and 
combustion rate. 

Generally, limited development in integrated engine predictive 
methods for direct injection engines has been due, in part, to individual 
component modeling for direct injection diesel fuelled engines that has 
been satisfactory for the purposes of modeling such engines. 



For example, the changing conditions inside the combustion chamber 
caused by combustion in a diesel engine has limited influence on the 
injection of diesel fuel due, largely, to the extremely high pressures at 
which diesel fuel is injected, so an integrated model that adjusts the fuel 
delivery to a cylinder in response to the cylinder pressure does not suffer 
to any significant extent by ignoring combustion of the fuel in the 
cylinder. That is, compared to the pressures established in the combustion 
chamber as a result of combustion events, the pressure at which diesel fuel 
is injected is many times higher and, therefore, less influenced by changes 
in combustion chamber pressure. 

However, ignoring this part of an engine system for a predictive 
model can be a problem when a gaseous fuel is used to operate high 
compression ratio direct injection engines. For diesel engines operating 
on gaseous fuels, the fuel injection rate depends in general on pressure in 
the combustion chamber; however the rate of change of cylinder pressure 
depends on the injection rate, as well as on the turbo-machinery and heat 
transfer components. Thus for a full system integration separate models 
for predicting injector behavior and combustion development are 
important. These models exist independently but have not been adequately 
integrated into a full system applicable to diesel engines, or engines with 
diesel engine compression ratios, where those engines operate on gaseous 
fuels. 

Most component models simulate aspects of an engine removed 
from the operation of the engine system as a whole. The interface between 
the modeled con^)onents and the rest of the system is based on engine test 
measurements or data. Therefore, the model is not a self-supported fiill- 
predictive model. The models are applied to understand the engine process 
for diagnostic purpose. 

As noted above, in injector computation, the model upstream 
conditions, such as rail pressure and temperature, are from measurements 
that, in general, do not consider fluctuations in pressure, which might be 
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quite significant in some cases. This is especially important when gas fuel 
injection is involved. In combustion models, the intake conditions are 
usually from test measurements on a functioning engine. There is little 
consideration of the interaction between the injection and the combustion. 
5 Likewise, in engine cycle simulation, the combustion process is often 
derived jfrom the heat release rate measured on an actual engine. This 
ignores the effect of the intake conditions on the combustion process and 
simplifies the interaction between injection, combustion and air intake and 
management processes. It is believed that the interaction between 

1 0 components in the engine system is very important and needs to be well 
understood in order to shorten engine development phase. However, it is 
also very complex. Most engine models do not reflect the con^lex 
interactions in engine system. Models are detached from the system 
reducing a model's predictive capability. 

15 Moreover, until recently computer power, has not been adequate to 

accurately integrate modeling components to predict overall engine 
behavior in a time frame useful for the purposes contemplated. 
Developments in computer processing power have provided the necessary 
computational capacity to model an integrated engine. 

2 0 Efforts to integrate sub-models together have, in general, not 

resulted in integrated fuel systems, combustion system and engine 
breathing systems (as well as emissions systems). In particular, modeling 
efforts have not been developed for such a predictive integrated system for 
a gaseous-fuelled engine. Also, while the focus of efforts to create engine 

25 predictive methods has been directed to the purpose of providing an 
engine diagnostics tool and predicting engine behavior, little effort has 
been made on the use of such methods for engine controls. 

The present method discloses an integrated engine prediction tool 
for controlling and predicting engine performance and emissions. 

30 
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Summary of the Invention 

The present invention provides for an integrated engine model and 
method for predicting and controlling a direct injection fuelled engine and 
5 especially a gaseous-fuelled direct injection engine. 

A further embodiment of the present invention includes A method 
is disclosed for controlling a direct injection intemal combustion engine 
by calculating output data from system data. The method comprises using 
system data which includes an injector geometry defined by a fuel injector 

10 in fluid communication with a combustion chamber of the engine, an 
engine geometry defined by the engine, an engine operating command 
based on a desired engine output, that is indicative of a commanded torque 
for the engine, an initial injector command, and an engine speed. The 
output data comprises an engine out value indicative of a brake torque 

1 5 delivered by the engine, and an output control parameter indicative of a 
cylinder pressure value. When the output data does not satisfy a 
predetermined relationship, it is recalculated with a determined subsequent 
injector command. This predetermined relationship compares the output 
data, the engine operating command and a demanded control parameter. 

20 The demanded control parameter is indicative of a maximum cylinder 
pressure value. When the output data satisfies the predetermined 
relationship, the injector is commanded to provide fiiel to the combustion 
chamber according to a last injector command. 

In a further embodiment of method, the system data can also 

2 5 comprise an ambient temperature, an ambient pressure and an engine 

speed, which in a preferred embodiment are measured prior to calculating 
the output data. In a further embodiment of the present method, the 
engine operating conmiand is converted to the commanded torque based 
on an operator selected pedal position. 

30 In a further embodiment of method, the demanded control 

parameter is determined from the oigine operating command. As well. 
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the output control parameter is indicative of calculated emissions and the 
demanded control parameter is indicative of a demanded emission. The 
calculated emissions and the demanded emissions are indicative of NOx 
and/or particulate matter concentration. 
5 A further embodiment of method discloses the predetermined 

relationship comprising comparing the output data with a previously 
calculated fuel consumption determined from the system data using the 
initial or subsequent injector commands. The output data further 
comprises the calculated fuel consimiption. To meet the predetermined 

10 relationship, the calculated fuel consumption is less than the previously 
calculated fuel consumption. 

In a further embodiment of method, the initial injector command 
comprises commanding initial start of injection, an initial pulse width and 
an initial rail pressure. The subsequent injector command comprises 

1 5 commanding a subsequent start of injection, a subsequent pulse width and 
subsequent rail pressure. And the last injector command comprises 
commanding a last start of injection, a last pulse width and a last rail 
pressure. 

A further embodiment of method discloses a variable geometry 
2 0 turbine. The system data may comprise variable turbine geometry 

positions. The injector commands include the last variable geometry 
turbine positions when the output data satisfies the predetermined 
relationship. When the output data does not satisfy the predetermined 
relationship, a subsequent variable geometry turbine position, with the 

2 5 subsequent injector command, is used to recalculating with the output 

data. 

In a further embodiment of the method the engine comprises an 
exhaust gas recirculation (EGR) system. The system data further 
comprises an initial effective EGR valve flow area, further comprising, 

3 0 when the output data does not satisfy the predetermined relationship, 

determining a subsequent effective EGR valve flow area with the 
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subsequent effective EGR valve flow area and the subsequent injector 
command. The output data being recalculated when the output data 
satisfies the predetermined relationship resulting in commanding the EGR 
system according to a last effective EGR valve flow area. 
5 In a further embodiment of the method, the output data is 

calculated using an injector module for calculating IM/CM data and 
IM/ECM data from an estimated initial cylinder pressure and IM system 
data. The IM system data is a subset ofthe system data. Also, a 
combustion module for calculating CM/ECM data from an estimated 

1 0 initial intake flow rate, the IM/CM data and CM system data is used. The 
CM system data is a subset of the system data. Also, an engine cycle 
module for calculating the output data from the CM/ECM data, the 
IM/ECM data and the ECM system data is used. The ECM system data is 
a subset ofthe system data. Prior to calculating the output data, a small 

15 loop between the combustion model and the engine cycle module is used 
to provide a converged intake flow rate. A big loop between the injection, 
combustion and engine cycle modules is used to provide a converged 
cylinder pressure, which then provides the cylinder pressure value. 

In a further embodiment of the present method, the IM/CM data 

2 0 can comprise an IM rate of injection and an IM start of injection. An IM 
fiiel flow and an IM start of injection can also be used. As well, the 
CM/ECM data can conqDrise a CM heat release rate and a CM start of 
combustion. 

The output data is further calculated using an emissions module for 

2 5 calculating emissions data from ECM/EM data, the IM/ECM data and the 

system data. The output control parameter is indicative of calculated 
emissions and the demanded control parameter is indicative of a 
demanded emission. 

Further embodiments of the present method include the IM system 

3 0 data comprising an ambient pressure, the injector geometry, the engine 

speed and the injector command. The CM system data can further 
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comprise the engine speed and the engine geometry. The ECM system 
data comprises an ambient pressure, an ambient temperature, the engine 
geometry and the engine speed. 

The method can be practiced wherein the engine is at least partially 
5 fuelled by a gaseous fuel and in a preferred embodiment the gaseous fuel 
is natural gas. The gaseous fuel can also comprise methane, ethanol and 
propane. In a further preferred embodiment the fuel is hydrogen. 

A method is disclosed for predicting output data generated by a 
direct injection internal combustion engine, comprising at least one of 

1 0 engine torque, fuel consumption, power and cylinder pressure. The 
method comprises selecting system data which comprises an injector 
geometry indicative of a selected fuel injector, an engine geometry 
indicative of the engine, an engine speed and an injector command. The 
method also comprises estimating initial data, comprising an initial 

15 cylinder pressure and an initial intake flow. The output data is calculated 
using the system data and a converged cylinder pressure and a converged 
intake flow. The initial data is used to provide the converged cylinder 
pressure from the initial cylinder pressure, and provide the converged 
intake flow from the initial intake flow. 

20 In a further embodiment, the system data can comprise either an 

ambient pressure and/or an ambient temperature while the output data can 
comprise NOx and/or particulate matter concentration. As well, the 
injector command may comprise a commanded rail pressure, a 
commanded start of injection and a commanded pulse width. 

25 In a further embodiment of the present method, system data can 

comprise a commanded variable geometry turbine position and an 
effective EGR valve flow area. When an effective EGR valve flow area is 
included, the system data further comprises an initial EGR flow. 

In a further embodiment the output data is calculated using an 

3 0 injector module for calculating IM/CM data and IM/ECM data from the 
initial cylinder pressure and IM system data. The IM system data is a 



subset of the system data. Also, a combustion module is used for 
calculating CM/ECM data from the initial intake flow, the IM/CM data 
and CM system data. The CM system data is a subset of the system data. 
Also an engine cycle module is used for calculating the output data from 
the CM/ECM data, the IM/ECM data and ECM system data. The ECM 
system data is a subset of the system data. With this, a small loop between 
the combustion model and the engine cycle module is used to converge to 
a converged intake flow. As well, a big loop between the injection, 
combustion and the engine cycle modules is used to converge to the 
converged cylinder pressure. 

Li a fiirfher embodiment of the method, the output data is fiirther 
calculated using an emissions module for calculating emissions data from 
ECM/EM data, the IM/ECM data and EM system data. The EM system 
data is a subset of the system data. The output data comprises the 
emissions data, which, in a preferred embodiment, is indicative of NOx 
concentration and/or a particulate matter concentration. 

In a further embodiment, the BCM/EM data comprises the 
converged cylinder pressure, the converged intake flow, an mtake valve 
closing time cylinder pressure and an intake valve closing time cylinder 
temperature. The intake valve closing time cylinder pressure and 
temperature are provided by the engine cycle module. 

In a fiirther embodiment of the subject invention, the IM/CM data 
comprises an IM rate of injection and an IM start of injection. The 
IM/ECM data comprises an IM fiiel flow and an IM start of injection. The 
CM/ECM data comprises a CM heat release rate and a CM start of 
combustion. The IM system data can comprise the injector geometry, the 
engine speed, the ambient pressure and the injector command. The CM 
system data comprises the engine speed and the engine geometry. The 
ECM system data comprises the engine geometry, the ambient pressure, 
the ambient temperature and the engine speed. 
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The method can be practiced wherein the engine is at least partially 
fuelled by a gaseous fuel and, in a preferred example, the gaseous fuel is 
natural gas. The gaseous fuel can also comprise methane, ethanol and 
propane. In a further preferred example the gaseous fuel is hydrogen. 

Brief Description of the Drawings 

FIG. 1 shows a flow diagram of an integrated engine model for use 
in engine control or predicting engine behavior. 

FIG. 2 shows a flow diagram demonstrating a small iteration loop. 
FIG. 3 shows a flow diagram demonstrating a big iteration loop. 
FIG. 4 shows the data flow path of the initial big loop as shown in 

FIG. 3. 

FIG. 5 shows the data flow path(s) for tfie subsequent big loops 
subsequent to the initial big loop of FIG. 4. 

FIG. 6 shows the data flow path wherein the integrated engine 
model is used for engine controls. 

Detailed Description of Preferred Embodiment(s) 

In the present disclosure, a method is taught to help with engine 
control, particularly for gaseous-fiielled direct injection internal 
combustion engines. The method predicts internal combustion engine 
behavior. 

Referring to Fig. 1, the flow diagram illustrates integrated engine 
method (lEM) 8. The flow diagram shows the integration of injector 
module (IM) 10, combustion module (CM) 12, engine cycle module 
(ECM) 14 and emissions module (EM) 16, each of which calculates output 
data based on provided input data where the output data is descriptive of 
the behavior of the injector, the selected combustion mode for ttie chosen 
fiiel, the engine behavior and the emissions generated by the engine of 
interest. System data 18 includes geometry parameters and operation 
parameters, to be described in greater detail below. System data 18 is 
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used to initiate operation of lEM 8. System data 18 lists all data used for 
lEM 8. Subsets of system data 18, IM system data 20b. CM system data 
21b, ECM system data 3 1 and EM system data 39, each provide system 
data 18 used for each module. 
5 As well, initial data, namely, initial cylinder pressure ( ) data 1 9 

and initial intake air flowrate (Q^ ). Exhaust gas recirculation flow rate 
{Qegr ) is also provided as initial data in those cases where exhaust gas 
recirculation (EGR) is used or considered. Q?„ and Q^^^ are, collectively, 
CM input data 21 where EGR is used and Qf„ alone where it is not used. 
10 Lines with arrowheads represent data paths for routing data 

throughout IBM 8 between modules 10, 12, 14 and 16, as well as routing 
system data 18, P^^ data 19 and CM input data 21 into lEM 8. Data paths 

20, 22 and 24b direct IM system data 20b, P^f data 19, P^, data 23, 

respectively, to IM 10. Data paths 26 and 28 direct IM/ECM/EM data 25 
15 and IM/CM data 27, respectively, from IM 10 as indicated. IM/ECM/EM 
data 25 is directed to ECM 14 and EM 16. IM/CM data 27 is directed to 
CM 12. IM system data 20b, a subset of system data 18 comprises, 
injector geometry data ( GEOi„j ), injector commands, in the embodiment 

shown rail pressure ( P^.^ ), commanded start of combustion ( SOI^ ), and 

2 0 commanded pulse width ( PW^ ), as well as and engine speed ( RPM^^ ) 

and ambient pressure of the engine of interest {P^^i)- 

Data path 24b originates from path 24 and directs data 23 to 

IM 10 after completion of an iteration of big loop as will be discussed 
below- 

2 5 Along with IM/CM data 27 and CM input data 2 1 (directed 

through path 40), as noted above, CM 12 also takes CM system data 21b 
through path 30c and generates CM/ECM data 29. After an iteration 
through ECM 14, as will be described below, CM 12 also uses ECM/CM 
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data 35 directed througji path 34 replacing CM input data 21 on 
subsequent iterations. 

ECM 14 uses IM/ECM data 25, CM/ECM data 29 and ECM 
system data 31 from paths 26b, 32 and 30a, respectfully, to provide 
data 23, ECM/EM data 33 and ECM/CM data 35 directed to IM 10, EM 
16, and CM 12 through paths 24b, 24a, and 34, respectfully. ECM system 
data 31 is a subset of system data 18. Path 36 provides engine output data 
38, which can be generated with each iteration or upon completion of a 
number of iterations as described below. Path 34 along with path 40, path 
28 and path 30c, as noted above, deliver ECM/CM data 35, CM input data 
21, IM/CM data 27 and CM system data 21b to CM 12. 

Path 24a, 26a and 30b direct ECM/EM data 33, IM/ECM/EM data 
25 and EM system data 39, respectively, to EM 16 from ECM 14, IM 10 
and system data 18, respectively, wherein EM 16 generates emissions data 
41 through line 42. EM system data 39 is a subset of system data 18. 

Along with IM system data 20b described above, other data, ECM 
system data 31, CM system data 21b and EM system data 39, are chosen 
from system data 18. ECM system data 31 includes engine geometry 
( (^^O^g ). ^Af„g , P^^ and ambient temporature ( ). EGR valve 
flow area ( A^^ ) and variable geometry turbine position ( FGr^„ ) may also 
be included where these features, EGR and variable geometry turbine 
(VGT) are used in the engine considered or controlled. Other variable 
may also be included depending on the engine, as would be understood by 
a person skilled in the art. 

Note while P^^^ and T^^ are components of system data 18 to 
both IM 10 and ECM 14, they can be assumed or calibrated values and 
factored into the computational algorithms in both IMIO and ECM 14 
thereby avoiding actually having to provide these as values. This could 
also be said of GEO„^ , GEO^ and RPM„^ , however, the flexftility of 
the system is limited by setting this parameters as constants. For the 
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purposes of this application, these three parameters will be considered 
provided as a subset the system data even if they are "hardwired" into the 
modules that use them. 

CM system data 21b and EM system data 39 include GEO^g and 

5 RP^feng ^ cases. Depending on such things as aftertreatment in the 

case of EM 16, and the adaptability of EM 16, other variables may be 
employed. 

As would be understood by a person skilled in the art, GEO^g is a 

reference defined by a number of variables including such things as the 
10 engine cylinder diameter, number of cylinders, piston stroke, connecting 
rod length, clearance height, compression ratio, strokes per cycle (two 
stroke cycle or four stroke cycle), number of intake and exhaust valves per 
cylinder, diameters of valves, valve lift and timing diagrams, number of 
compressors and turbines (if any — although, its assimied for the purposes 
15 of this disclosure that these are feature of the engine), size and properties 
of compressor, size and properties of any turbine. 

As would be understood by a person skilled in the art, GEOf„j is a 

reference defined by variables including sizes of injector needles, sizes of 
plungers, sizes of flow channels and reservoirs, sizes of injector holes and 

2 0 number of inj ector holes. 

Output data, engine output data 38 and emissions data 41, includes 
power (Pfr^) and torque {Tf^rake)* economy (JBiSFC), exhaust 
emissions ( NO^ , PM ), intake air flow ( Q^^ ), P^i and temperature and 
pressure of the turbine inlet ( 7],. , P,, ) and the compressor outlet ( T^^ , P^^ ), 

2 5 turbine speed ( v,^ ). Where used, EGR rate ( ) can also be provided 
with engine output data 38. As would be understood by a person skilled in 
the art, variables such as exhaust flow {Q^ ), where important can be 
found. Heat release rate ( HRR ), an intermediate parameter out of ECM 
14 originally delivered to ECM with CM/ECM data 29 from CM 12, can 
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also be used by the engine cycle model to generate engine performance 
parameters, such as P,^^ and r,„^, BSFC, Q^.Q^, P^, 7;.^ and 

TMICM data 27, IM/ECM/EM data 25 and CM/ECM data 29, 
includes rate of injection (ROI\ determined start of injection (50/ ), in 
5 the case of IM/CM data 27. In the case of IM/ECM/EM data 25, SOI and 
rate of fuelling are mcluded. In the case of CM/ECM data 29, 

HRR and start ofcombustion (500)15 provided. These data sets, 27, 25 
and 29, form the basis of two iteration loops used to converge on a P^^ 
data 23 set and a ECM/CM data 35 set that more accurately reflects and 

10 refines initial input data, P^^ data 19 and CM mput data 12. ECM/CM 

data 35 includes the current iteration of Q.^ , referred to as , and may or 
may not include the current iteration of Q^^ , referred to as Q^^ , 
depending on whether this is a feature of the engine used. 

ECM/EM data 33, includes intake valve closing time cylinder 

1 5 pressure ( P^^^ ), intake valve closing time cylinder temperature ( T^^^ ) as 
well as Q^^ , P^^ , and, where a feature of the system, Q^^^ . With the flow 
rates of fiiel, air, and EGR, pressure and temperature at the intake valve 
closure, lEM 8 module, EM 16 calculates emissions data 41, such as NOx 
emission. 

2 0 Referring to figs. 2 and 3, small loop 98 and big loop 1 1 8 are 

iteration algorithms which are used to establish initial input data of 
sufficient accuracy to allow lEM 8 to begin self supported operation and 
generate engine output data 38 and emission data 41, as will be discussed 
below. 

2 5 Referring to fig. 2, small loop 98 iteration, initial data, CM input data 2 1 , 
CM system data 21b and IM/CM data 27, are directed through flow paths 
40, 30c, 28, respectfiilly, to CM 12. CM 12 calculates CM/ECM data 29 
through path 32 to ECM 14 where, using ECM system data 31, ECM/CM 
data 35 is provided and directed back to CM 12 in a closed start-up 
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algorithm. ECM/CM data 35 for each iteration, 0,^ (and Q^^ where EGR 
is used), is compared with the previous value provided for ECM/CM data 
35, Q"~* (and 0^;' where EGR is used), as seen in small loop 
convergence step 200. 
5 In fig. 3, big loop 1 18 iteration shows the interaction between big 

loop 118 and small loop 98. Big loop 1 1 8 includes IM system data 20b 
and F^i data 19 and, eventually, the current iteration value of from 

data 23, namely, P^ , all of which are directed to IM 10 through paths 
20, 22, 24b. Resulting IM/CM data 27 is delivered through path 28 to CM 

10 12. CM 12 uses this data, CM system data 21b, CM input data 21 and, 
eventually, ECM/CM data 35, delivered through paths 30c, 40, 34 to 
generate CM/ECM data 29. CM/ECM data 29 is delivered through path 
32 to ECM 14. ECM 14 performs iterations of small loop 98 to generate 
Q"„ (and, where appropriate g^J until this value is converged as 

15 discussed below after which, ECM 14 generates P^ data 23 delivered 

through path 24b to IM 10, after conqjleting big loop convergence step 
202. 

The data flow paths for big loop 1 18 are shown in fi^. 4 and 5. 
Referring to fig. 4, initial flow path 160 is provided wherein big loop 1 1 8 

2 0 starts from IM 10, goes through a gas flow loop, returning to IM 1 0. 

Subsequent big loop flow path 162 is shown in fig. 5. Here, IM 10, CM 
12 and ECM 14 are cycled through sequentially. 

In practice, referring to fig. 1, each of the modules 10, 12, 14 and 
16, requires input data to generate results indicative of that aspect of the 

2 5 engine under consideration. The input data for individual modules is 

dependent on such things as the control accuracy sought, the demands of 
the user or operator and the time available for open loop control or engine 
behavior prediction. In general, certain data is required in order for the 
controls and engine prediction to be practiced, as would be understood by 
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a person skilled in the art. It is possible to supplement the controls and 
engine prediction algorithm with additional data depending on the needs 
of the user, the engine in question and the accuracy sought regarding the 
resulting data. 

5 An initial estimation is required to begin the subject method. The 

disclosure utilizes two iteration processes demonstrated in figs. 2 through 
5, each of which limits the impact of this initial estimation by limiting the 
estimation to data that is relatively easy to initially estimate and using big 
loop 118 and small loop 98 to correct for inaccuracies in the initial 

1 0 estimation by performing a converging step on this estimated values. 

JEM 8 uses system data 1 8, namely, geometry data, GEO^^j , 
G^O^,, operational data, SOI^,, PW^, RPM^^, T^, P^^, 
along with initial data P^, data 19 and CM input data 21 to generate 
output data, engine output data 38 and, if of interest, emissions data 41. 

15 CM input data 21 may be limited to alone where no EGR is used. 

Also, where used in an engine control, the subject method may not, per je, 
use GEO^„j , GEO^^ , but calculate output data in light of GEO^^j , 
GEO^^^ . That is, these geometry 'Variables" may not be variables in an 
engine controller as such an engine would maintain the same GEO^^j and 

20 GEO^„g throughout. However, flexibility is also available to allow for 
variable valve actuation resulting in a variable GEO^^j and variable 
compression ratio resulting in GEO^„^ . Each would be a difficult variable 
to map allowing lEM 8 to provide real time mapping for engine controls. 
'^omb » f'amb eng ^an bc mcasurcd and sent to a controller 

2 5 for control of an engine. These variables could be predetermined for 

engine behaviour prediction. Operational data, P^^,, SOI^, PW^.zre 
operator conunanded variables indicative of the demands required by the 
operator for the engine as translated by a controller or predetermined in 
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view of a sought engine behaviour prediction. 

As would be understood by a person skilled in the art. additional 
variables can be added depending on the engine controlled or the engine to 
be predicted. Exan5>les include KG7),„ and as shown in Fig. 1, 

5 reflective of an engine with optional components^ VGT and EGR. Other 
variable would need to be included where used with the engine in 
question. 

lEM 8 relies on generation of IM/ECM/EM data 25, CM/ECM 
data 29, ECM/EM data 33 by models IM 10, CM 12 and ECM 14 to 

10 provide engine output data 38 and emissions data 41 . Each one of these 
data sets, 25, 29 and 33, is dependent, however, on P^, data 23, and 
ECM/CM data 35, generated by iterations performed with big loop 118 
and small loop 98. That is, IM 10 requires converged values for data 
23 on subsequent cycles of big loop 1 18 in order to generate a data set for 

15 useful to provide the output data, 38 and 41. Likewise, CM reqmres 

converged values for ECM/CM data 35 on subsequent cycles of small 
loop 98 in order to generate a data set for Q^, (and Q^^ where EGR is 
used) that is useful both for deriving ou^ut data, 38 and 41 , but also for 
supporting convergence of i^, data 23 through big loop 118, 

2 0 ECM 14 and EM 16 need initial data in the method taught that is 

not available through the start-up of the method. Before ECM 14 can 
generate any useful data, it requires CM/ECM data 29. Likewise, EM 16 
requires ECM/EM data 33 prior to generating emissions data 41. As such, 
there are steps within the taught method that must be performed before the 

2 5 integrated system will perform as a whole. Initially, IM 1 0 and CM 1 2 
operate with ECM in big loop 1 1 8 and small loop 98 within lEM 8, as 
described below, to generate the first sets of generated data, data 23, 
CM/ECM data 29. ECM/EM data 33. ECM/CM data 35. 

Referring to figs. 2 and 3, big loop 118 and small loop 98 initiate 
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lEM 8 by stepping through two loops. The method first performs 
iterations through small loop 98 convergmg on a (and Q^^ where this 
is used). Once, completed, big loop iteration is continued, converging on 

5 Initiating big loop, an estimate of P^, is made, namely P^, . This 

allows IM/CM data 27 to be generated. CM 12 is then able to take 
initially available estimated CM input data 21, (and Q^^ where EGR 

is used) and CM system data 21b along with IM/CM data 27 to initiate 
small loop 98. CM 12 generates CM/ECM data 29 which in turn allows 
1 0 ECM to use this data plus ECM system data 3 1 to provide a first iterated 
value of , specifically, Ql , and for each subsequent iteration Q^^ . 

Again, where EGR is used, a first iterated value is provided, namely, 

Qlg^ , and for each subsequent iteration Q^^^ . 

Convergence step 200 is performed wherein the relationship shown 
15 is considered. If true, then the convergence step is complete. More 

generally, this step involves a comparison of the first iteration value of 
Qin » Cm > to . In general however, the iteration steps involve the 
comparison of current Q„ , fi^ , to the previous value of , Q"~^ . The 
absolute value of the difference between Ql^ and Q"^'^ , AQ.^ , is compared 
20 to an error, , predetermined as acceptable for the application. An 
additional iteration is performed when 

|Aa„|=|a:-er'|>^e i 

and Qi„ is set such that = Q"^ when 

|AGj=|e;-er|^^e 2 

25 At this point, small loop iteration is conq)leted for the first cycle of big 
loop 118. 

For the purposes of this application Q^^ is the converged value of 
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defined as when \Q^„ -0-^1 ^ Sq. 

These values are directed back through CM 12 for subsequent 
iteration when \q^„ - Q"~' \> Sq, CM 12 again calculates new CM/ECM 

data 29, and 50C from IM/CM data 27, 50/ and ROI , and 

5 GEO^^ and RPM^^ , all of which remain the same through each iteration 
of small loop 98. 

The same process is performed to determine Q^^ where EGR is 

used. That is, the absolute value of the difference between O" and O""* , 
^Qegr » is compared to an error, Sq^^, predetermined as acceptable for the 
1 0 application. An additional iteration is performed when 

and ^ is set such that Q^^ = Q^^ when 

For the purposes of this application is the converged value of 
15 defined as Q"^ when |e;, < e^^, . 

Upon completion of this step where Q.^ = Q"^ (and Q^^^ = Q^^ , 

where appropriate), with the initial iteration cycle through small loop 98 
complete, iteration through big iteration loop continues. Here the resulting 
first iteration of P^, , in this case, but generally, P^ , is sent througji 

2 0 convergence step 202 wherein the relationship shown is considered. If 

tme, then the convergence step is complete. More generally, this involves 
a comparison of the first iteration value of P^^ , P^^ , to P^^ . In general 
however, the iteration steps involve the comparison of current i^,, to 
the previous value of P^ , P^;^ . The absolute value of the difference 

25 between P^^ and P^;^ , AP^, , is compared to an error, s^, predetermined 
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as acceptable for the application. An additional iteration is performed 
when 

K./|=kc;/-^;r|>^p 5 

and is set such that P^^ = P^^ when 

At this point, iteration of big loop 118 is completed where and Q„ 
(and Q^^^ where EGR is used) are determined allow for completion of 
lEM 8. Note that each iteration of P^, where 




7 



10 requires a further set of iteration steps through small loop 98 as described 
above generating an increasingly accurate value for Q^^ . 

For the purposes of this application P^^ is the converged value of 
defined as P^ when |p;, -/>;r^| < s,. 

Note also, that convergence considering |aP^;| and \M2in\ 
15 determined by equivalent means of comparing iterative values of both i^, 
and Q^^ resulting in a determination that the value has converged. These 
include a direct comparison of the iterative values without determining the 
absolute value or the difference in the values. As would be understood by 
a person skilled in the art other methods are also available for assessing 
2 0 convergence of a sought value. 

Examples of suitable values for Sq , €q^^ , Sp is less than 2% of 

the resulting Q"^ , Q"^^ , P^, . The value can be much higher or much lower 
depending on the needs of the operator. 

The number of iteration through either big loop 1 1 8 or small loop 
25 98 is, in general, determined by the user and the time available for each 
iteration as balanced against the accuracy sought for the result for the 
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application contemplated. This is provided by the value predeteirnined for 
each of Sp and Sq , in example method provided. 

An advantage of the ta\ight method is that the initial estimated data 
is limited to P^^^ and (plus, where EGR is used, ). As these two 

5 values are relatively easy to estimate for engine controls or engine 

behaviour prediction, as would be understood by a person skilled in the 
art, the efficacy of the control or accuracy of the prediction is determined 
in large part to these estimations and the accuracy of each module. 

While generally a function of time (typically expressed in degrees 
10 of crank angle rotation), may also be estimated as a single constant 

pressure, generally resulting in accurate results after completion of a few 
iterations of big loop 1 18 through subsequent big loop flow path 162. In 
the case of engine control, P^^ could come from the last completed 
predicted data set. 

15 Considering each module, IM 10 uses, along with P^^ data 19 in a 

first iteration and P^ in subsequent iterations, some of system data 18, 

namely, GEO,,^ , , PW^, , RPM^^ , and P^^, to generate 

data 27, namely ROl and SOI . Each one of IM system data 20b values 
is provided by the user as a commanded or desired or predetermined 
2 0 value. 

Turning to ECM 14, however, as well as a selection of data from 
system data 18, namely, GEO^^ , T^^^ , P^^^ and M*M^^ , ECM 14 also 
needs CM/ECM data 29, HRR and SOC, and IM/ECM/EM data 25, 
SOI and . These are relatively difficult values to estimate initially as 
25 they presume values dictated by an engine's upstream behaviour, i.e., 
actual SOI and not SOI^ , by way of example. However, CM 12 can 
provide such data along with IM 10. CM 12 uses CM input data 21, (and 
QLr where used). Along with this value, CM 12 uses IM/CM data 27 as 
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well, namely, SOI and ROI , and system data 18, namely, GEO^^ and 

RPM^^ , to provide a first iteration of HRR and SOC . Tliis commences 

small loop 98 wherein, HRR and SOC provided to initiate ECM 14 as 
described above. 

5 Referring again to fig. 1 together with figs. 4 and 5, data flow 

through initial big loop flow path 160 and subsequent big loop flow path 
162, see fig. 5, illustrate the derivation of ECM/EM data 33 used by EM 
16 along with EM system data 39 IM 10 generated IM/ECM/EM data 25, 
Qj^^^ and SOI , to generate emissions data 41. As none of IM 10, CM 12 
10 and ECM 14 use EM 16 generated data, EM 16 as a module can be 

removed from the iteration process and brought back in when P^^ and 

(and Q^^^ where EGR is used) are determined. 

The result of lEM 8 is the generation of emissions data 41 and 
engine output data 38, namely, NO^, PM , P,^, r,,,,,, BSFC, ft,, P^,, 

15 v,„^ , T^^ .T^n Pco^ Pti and where EGR is used. Other data can be 
generated depending on the engine considered and the ECM, CM, IM or 
EM modules used. The resulting values can be compared to the expected 
values and a third iteration process can be performed converging such 
things as NO^ , PM , P,,„,,, r,,„j^, BSFC choosing different IM system 

2 0 data 20b and different ECM system data 3 1 until the result sought is 

achieved. In effect, this is the iteration process used in practice where an 
operator will continue to demand, by way of example, more power 
contmuing to adjust values such as P^^^ , SOI^^^ and PW^ where these 
variables can be adjusted. For example, a driver commanding ever 

2 5 increasing power from the engine will adjust the pedal position at a rate 
sought to meet his or her needs. As would be understood by a person 
skilled in the art, the throttle translates into a given SOI^ and PW^ , in 
effect, and may in some engines allow for a change in i^,.; or 0 ^ , where 
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used. These changing variables allow the engine to drive to the demands 
sought of the driver in this example. Similarly, a sought emissions level 
can be met by considering the results provided by EM 16 after processing 
a first set of P^.^ , SOI^ and PW^ and then adjusting these to meet the 
5 emissions sought while still taking into consideration the other demands of 
the operator. In the case of controller, once the expected output data 
matches the actual output data with ECM output data 38 and emissions 
data 41, the controller will command the engine to operate with the P^^.^ , 
SOI^ and PW^ determined. Further discussion is provided below. 

10 The engine prediction application will provide output data 

indicative of the engine behavior when that data meets the goals of the 
JEM operator in some respect. 

The following table provides an outline of the types of variables 
that can be considered to provide engine output data 38 and emissions data 

15 41 for a particular engine design: 



Name 


Impact on: 


Injector needle diameter 


SOI , ROI , Pulse Width (PW), Q^^^, 


Injector tip hole size/number of 
holes 




Injector control chamber plunger 
diameter 


SOI,ROI,PW,Q^,,, 


Injector needle/plunger travel 
length 


ROI,PW,Q^^^, 


Injector needle/plunger springs 

mass/stiffhess 


SOI,ROI,PW,Q^,, 


Injector sac volume 


SOI.PW 


Combustion chamber clearance 
height - compression ratio 


P^.BSFC 


Turbo-charger turbine size 


Qin > ^turb jintake manifold pressure 

( ^im )» ^cyt > Qegr 


Turbo-charger compressor size 


Qin * ^turb > ^im » Qegr 


EGR configuration 


Q^,BSFC,NO„PM 


EGR valve size 


Qegr 
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iiilajvw CAllaUbl. ValVcS UXIUIlgy 11X1^ 

effective area 


Effective Volvime ( ), Q.^ , BSFC 


After treatment devices - pipeline 
layout 


Exhaust temperature( ), 


Command SOI & PW 


SOI ,ROI,PW, Q^^, , NO, , PM 


Fuel rail pressure 


SOI , ROI , PW , , NO, , PM 


EGR valve opening 


Q^„NO,,PM 


EGR coolant temperature 


EGR temoeraturer T \ NO PM 


EGR cooler size 




Intake/exhaust restrictions 


Flow restriction of intake and exhaust 
( Ar ), Q^^ , air to fiiel ratio 

(X\BSFC NO PM 


Elevation/environment temperature 


P*™*. . 7*™*. A, BSFC , NO, , PM 


VGT 


tsr,Q^,X, BSFC , NO, , PM 


Waste gate opening 


Ar , A. BSFC , NO, , PM 


Valve openings of after treatment 
devices 


Distribution of exhaust flow and 
temperature 



Note that for gaseous-fuelled engines, a pilot injection feature may 
be in^lemented within IM 10 with CM 12 for managing natural gas / pilot 
fuel combustion typical of a type of gaseous-fuelled direct injection high 
5 compression ratio engine. The resulting CM/ECM data 29 can be 

generated to support the remainder of the method as would be understood 
by a person skilled in the art. 

Referring to Fig, 6, the subject invention is described with one 
embodiment for use in an engine control environment. Here, operator 
10 command 300, from an operator or driver, provides an engine operating 
command. In this example, pedal position, a typical from a driver, is 
translated by a look-up table, in ECU 302, into torque command, . 
™s is directed with engine speed, RPM^^ , to ECU 302 from 
engme 305, see data file 304, to Supervisory Control Unit (SCU) 308, 
15 collectively referred to as SCU input data 306. SCU 308 provides, from 
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^cmd » operating boundary data 3 10 for engine operation to ECU 302. 
These include maximum cylinder pressure (P™ ), and can include 
demanded particulate emissions ( PM^ ), demanded NOx emissions 
( NO^"^ ), demanded temperature at turbine inlet ( T/'^ ) and demanded 
5 temperature at compressor outlet ( ). A modified torque can also be 

included where torque is to be kept within a certain range. 

System data is then provided to lEM 8 to allow it to calculate 
output data. Specifically, ECU 302 provides to lEM 8 a subset of input 
data 312, namely an initial injector command, in this case, an initial start 
10 of injection ( SOIl^ \ an initial rail pressure ( P^^,., ), an initial pulse width 
{PW^\ and and RPM^^ firom SCU input data 306. Optionally, an 
initial variable geometry turbine position ( YGT^^ ) and initial effective 
EGR valve flow area ( A^^^) can be provided where these variables reflect 
features of the subject engine. lEM 8 is also fed or stores base data 314, 
15 namely, GEO,^ , GEO^^ , T^^ and P^^^, the latter two variables provided 
by appropriate sensors. Its also possible the engine would be calibrated 
for a given T^^^ and , in which case, an actual measurement or 
estimate of T^^^ and would not be needed. 

With system data, namely input data 312, and base data 314, lEM 
20 8 can then provide output (predicted) data 316. For the purposes of 

control, output data should generally include engine out data including 
brake > ^ output control parameter which includes i^; and can include T^^ , 

> , v,„^ and PM and can include fiiel consumption, usually 
BSFC. These values are then compared to sought torque, r^^ , and any 
2 5 operating boundary data dictated by SCU 308 in data 310. In general, the 
demanded control parameter, CONP"^ , which is a fimction of some 
operating boundary data 3 10, is compared to the generated control 
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parameter provided in a first iteration, CONP'* . Two outcomes result 
through step 318. Where the relationship sought between CONP^'^ , 
provided through ECU 302 as parameter data 319, and CONP" and 
and T^^^^^ does not satisfy the conditions provided for in step 318, namely: 

CONP" <. CONP""^ 9 

where 

CONP" =f(NO,,PM,P^,T^.T,) 10 
CONP'''^ = f(NOi'^ ,PM'^''\P^\T^^\Tf,"^) 11 
10 while the determined fuel consumption, BSFC , satisfies a 

minimum value to meet the conditions above ( min (55FC)), a further 
iteration is performed through step 320 to determine a new predicted data 
set 316. Specifically, a convergence algorithm for deriving a multi- 
dimensional solution, such as simplex 322, is used to generate a new set of 
15 input data 324 which is, along with and RPM^^^ , delivered to lEM 8 

to allow it to generate predicted data set 316. Each n iteration through 
input data 324 provides a subsequent injector command, again, in this 
case, start of injection {SOI^ , n= 1,2,..., last), rail pressure (P^,, , 

n=l,2,...,last), pulse width (PW^ , n= 1,2,..., last) and variable geometry 

2 0 turbine position ( VGT^^ ,n=l ,2,. . .,last) and effective EGR valve flow area 

( A^g^ ,n= 1 ,2, . . . ,last), where appropriate. 

Note, in general, CONP'''^ is a function of P^"" as the cylinder 

pressure would always be constrained by the physical limitations of the 
cylinder and need to be monitored. The other variables noted in equation 
25 11 can be added depending on the engine used and the variables of 
interest. Also, while the SCU 308 provides the parameters used to 
determine CONP"^"^ , SCU 308 can also provide the relationship directly 
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to ECU 302 or lEM 8 for use in step 318. 

Eventually, step 318 will be satisfied resulting in step 326 through 
to ECU 302, where the last values generated for input data 324 are 
delivered to allow ECU 302 to then direct engine 305 to operate using the 
5 last derived injector command, namely, a start of injection ( SOI^^ ), a rail 

pressure ( P^f) and a pulse width ( PW^^^) and a variable geometry turbine 
position ( VGT^ ) and an effective EGR valve flow area ( ), where 
appropriate, where each is set to the respective last value determined for 
SOr^, PZu, PW^ and VGT'^ and where appropriate. That is, 

10 SOI^^SOI'Z 12 

PW^=PW^ 13 

P^=P':^, 14 

VGr^^VGTX 15 

= A'^ 16 

15 In practice, such variables as the rate of change of pedal position 

are taken into consideration. The iterative process discussed above can be 
performed over a series of different pedal positions, which should provide 
a rate of change in the pedal positions being commanded of the operator. 
The comparison made between the operator commands or the 

2 0 needs of the operator as compared to the conditions determined by lEM 8 
as a possible means of meeting those needs is approximated within a range 
available for engine operation. That is, while the operator generally 
requires a given torque for his or her purposes, limitation of the engine 
physically or limitation set by SCU 308 for other reasons may prohibit the 

25 sought torque under a set of engine commands. As such, ECU 302 may 
use, as noted above, additional boundary data such as a modified torque 
able to approximate the operator demands and still operate within the 
physical boundaries of the engine and vehicle or turbine design. 
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In practice, ECU 302 can include all the components physically 
independent in fig. 6. In other words, ECU 302, SCU 308, lEM 8 and 
Simplex 322 can be physically separated as suggested, but would likely be 
separate routines within the ECU 302 itself. 
5 A further embodiment of the control strategy can include 

diagnostic capabilities that use an engine condition value indicative of a 
condition in the engine to monitor actual sensors in the engine. Here 
variations between predicted operating conditions and measured operating 
conditions can identify failures in such things as emissions sensors, EGR 

10 flow monitors, knock sensors and other combustion quality sensors, all of 
which could be fed, for such diagnostic analysis, to ECU 302 with 
predicted data 316 and input data 324 satisfying step 318. Faulty sensors 
can then be discovered and ignored by the ECU 302 for future cycles. On 
board diagnostic functions can be alerted as to the existence of faulty 

15 sensors and appropriate actions taken according to on-board-diagnostic 
regulations. As well, replacing those faulty sensor inputs with the IBM 8 
generated values can be xised for limited model-based engine control. 

Still a further embodiment provides for a control strategy that 
works with in-cylinder pressure sensors. Comparisons can be made 

2 0 between the predicted values of cylinder conditions at or near start of 
combustion as predicted by lEM 8 and measured in-cylinder pressure 
sensors. This can help supplement controls of variable used to adjust start 
of combustion or combustion phasing. 

Similarly, perturbations in variable rate injectors can be predicted 

2 5 and controlled by lEM 8 helping to enhance the dependability of such 
injectors which can help to capitalize on the value of such injectors. 

For controls or predictive models, as would be understood by a 
person skilled in the art, the variables used by JEM 8 initially such as input 
data 312 or IM system data 20b can be other data that effectively provides 

30 (he same information. For example while the injector commands of input 
data 312 include SOI^.P;^ , and PW^^ , the injector commands can 
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also be reflected by an initial SOI^ , ROI^ , and Q^^ with the same or 
similar resulte. Other variables can be used to determine and control fuel 
injection and combustion characteristics. 

5 Example of lEM Operation 

lEM 8 integrates each of IBM's 4 modules, IM 10, CM 12, ECM 
14 and EM 16, consisting of a number of data files, which can be held in 
independent directories depending on their properties and functions. 
10 An exan^le of the subject method for engine behaviour prediction 

is provided. The example does not limit ttie scope of the invention but is 
meant only for the purposes of illustrating an embodiment of the 
invention. For IM 10, the data files used include control files, P^^ and 

files, command files, and output files. The control file defines SOI^^ 
15 and PW^ of injection signals and are, therefore, a subset of IM system 
data 20b. By way of example, the following is typical of this file 
structure: 



12 



0.0 


0.0 


0.0 


0.0 


1.000 


0.0 


0.0 


0.0 


1.010 


1.0 


0.0 


50.0 


1.100 


1.0 


0.0 


50.0 


1.110 


1.0 


0.0 


0.0 


1.140 


1.0 


0.0 


0.0 


1.150 


1.0 


0.0 


9.0 


1.700 


1.0 


0.0 


9.0 


1-701 


1.0 


2.0 


9.0 


3.000 


1.0 


2.0 


9.0 


3.010 


0.0 


0.0 


0.0 


10.000 


0.0 


0.0 


0.0 



2 0 which defines the time (first column in ms) and commanded opening and 
closing timings (second column) where zero means closing, one means 
opening. In the example, referring to IM system data 20b, P^^ is a 
constant, an example of the value for a gaseous fuelled engine being 
185e5 Pa. 
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Additionally IM 10 requires, P^i data files, files that hold data 
19. By way of example, a backpressure file could look like: 

1440 0 . S 

0.0 41.173 745.985 

0.5 42.314 751.362 

1.0 43.493 756.802 

1.5 44.711 762.306 

2.0 45,969 767.871 

2.5 47.269 773.498 

3.0 48.611 779.186 



717.0 35.070 715.074 

717.5 36.005 720.065 

718.0 36.971 725.121 

718.5 37.970 730.240 

719.0 39.002 735.424 

719.5 40.069 740.671 



5 The first line provides total number of P^^ data points in the file 

and the crank angle interval between data points. From the second line 
forward, as an initial reading, the crank angle (CA) is provided in the first 
column, and P^^ , measured in bar, is provided in the second column. In 
this example, P^f starts from zero degrees C A, which corresponds to zero 
10 timing of the control file. A reading of 0 CA of the F^^, file does not mean 
TDC of an engine cycle, but means the start point of the control file. The 
injection timing corresponding to an engine cycle can be varied by 
changing the starting point of P^, . 

Note as well, that P^^ is, as mentioned above, an estimate of the 

15 backpressure in the system at issue. In general, this estimate can come 
from data generated on an actual engine with similar operational and 
physical specifications. Also, could be set equal to the cylinder 
pressure determined from a previous converged pressure value (as is the 
case for , see below). However, where such a file is not available, a 

2 0 constant pressure, while a crude estimate, can be used for P^j . After 

iterations of big loop 118, P. will provide a more indicative measure of 
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backpressure during the course of an engine cycle on the engine at issue. 
An example of a constant P^^ could be: lO.OeS kPa. 

The command file for IM 10, IM system data 20b, defines the 
injector settings and operational conditions. These are operator driven and 
5 provide parameters of the system considered or system controlled. The 
data includes the variables and constants found in IM system data 20b and 
IM input data 19 as well as typical specifications required for the modules 
to operate and defining the boundaries of the iteration, for example, how 
long to run the model, the maximum time step allowed, where to output 
10 results. 

Also, measiared or predefined values for RPM^^ and measured 
values for P^^ are used. An example of RPM^^ is 1740 and is 
1 .0e5 Pa. This will allow big loop 1 1 8 to commence by initiating IM 10. 
Also, GEO^^j file could be referred to where the file would 
15 include, typically, information such as: 



Needle diameter 


[mm] 


Needle mass 


[gram] 


Needle stiffness 


[N/m] 


Plunger diameter 


[mm] 


Plunger mass 


[gram] 


Plunder stifftiess 


[N/m] 


Injector hole effective flow area 


[mm^2] 


Injector hole L/D ratio 




Needle seat angle 


[degree] 


Needle spring mass 


[gram] 


Needle spring stiffiiess 


[N/m] 


Plunder spring mass 


[gram] 
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Typical units are provided in the right hand column. 

RPM^^phxs GEO,„j, all 

noted above, typical output data files can be generated firom IM 10. The 
main file to be used by lEM 8 is the formatted file, an example of which is 
provided below: 



1740 . 0 

-30 .0 

-29.0 

-28.0 

-27.0 

-26.0 

-25.0 

-5.0 
-4 .0 
-3.0 
-2.0 
-1.0 
0.0 
1-0 
2.0 
3.0 
4.0 
5.0 



-4,8 16 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOEi-OO 
O.OOOOOE+OO 
O.OOOOOE-l-00 



O.OOOOOE-i-OO 
0-OOOOOE+OO 
O.OOOOOE-i-OO 
O.OOOOOE-kOO 
O.OOOOOE+OO 



0 .OOOOOE+OO 
O.OOOOOE+00 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 



0. OOOOOE+OO 
0. OOOOOE+OO 
O.OOOOOE+00 
O.OOOOOE+00 
0. OOOOOE+OO 



0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 0. OOOOOE+OO 



0. OOOOOE+OO 
0.26492E-06 
0,18489E-05 
0.38032E-05 
0.50160E-05 
0.56768E-05 
0.61400E-05 
0.65121E-05 
0.68056E-OS 
0.68804E-05 
0.6901BE-05 



0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0 .OOOOOEi^OO 
0 .OOOOOE+OO 
O.OOOOOE+00 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 



O.OOOOOE+00 
O.OOOOOE+00 
0. OOOOOE+OO 
0 .OOOOOE+OO 
0 .OOOOOE+OO 
0. OOOOOE+OO 
O.OOOOOE+00 
0. OOOOOE+OO 
0. OOOOOE+OO 
O.OOOOOE+00 
0. OOOOOE+OO 



0. OOOOOE+OO 
0. OOOOOE+OO 
O.OQOQOE+00 
0 . OOOOOE+OO 
0 , OOOOOE+OO 
O.OOOOOE+00 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
O.OOOOOE+00 
0. OOOOOE+OO 



65 


0 


0 


OOOOOE+OO 


0 .OOOOOE+OO 


0 


OOOOOE+OO 


0 


OOOOOE+OO 


66 


0 


0 


OOOOOE+OO 


0 .OOOOOE + OO 


0 


OOOOOE+OO 


0 


OOOOOE+OO 


67 


0 


0 


OOOOOE+OO 


0, OOOOOE+OO 


0 


OOOOOE+OO 


0 


OOOOOE-fOO 


68 


0 


0 


OOOOOE+OO 


0 .OOOOOE+OO 


0 


OOOOOE+OO 


0 


OOOOOE+OO 


69 


0 


0 


.OOOOOE+OO 


O.OOOOOE+00 


0 


OOOOOE+OO 


0 


OOOOOB+00 


70 


.0 


0 


.OOOOOE+OO 


0. OOOOOE+OO 


0 


.OOOOOE+OO 


0 


OOOOOE+OO 



The first number being RPM^^ , followed by SOI (measured in degree 

10 of crank angle rotation relative to top dead center), and total fuelling 

(measured in mg/shot). The first column is degrees of crank angle rotation 
and the second column is ROI (measured in [kg/ca] in this case). is 

derived fi-om an integration of ROI data as would be understood by a 
person skilled in the art. The gap in data between crank angles —25 and -5 
15 and 5 and 65 is given for economy. The data file would not include these 
gaps. 

The result is IM 10 generating IM/CM data 27 and IM/ECM/EM 
data 25. 
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Having generated IM/CM data 27, CM 12 uses this data plus its 
own data file, CM input data 21, (and Q^, if EGR is a consideration in the 
engine considered) as well as CM system data 21b» RPM^^ and GEO^^ . 
An example of initial CM input data 21 with EGR is = 1500.0 kg/hr 
5 and Q^g, = 300.0 kg/hr for natural gas. 

CM 12 provides CM/ECM data 29 in a file for use by ECM 14. 
An example of such a file is: 



•20 


0 


-19 


0 


-18 


0 


-17 


0 


•16 


0 


-15 


0 






-5 


5.36444E-08 


-4 


1.42945E-05 


-3 


0.00012872 


-2 


0.00034991 


-1 


0.000618775 


0 


0.000939128 


1 


0.001366183 


2 


0.00197213 


3 


0.002863935 


4 


0.00416602 


S 


0.006229198 






95 


0 


96 


0 


97 


0 


98 


0 


99 


0 


100 


0 



1 0 The CA rotation relative to top dead center is provided in the left 

hand colunan and HRR is provided in the right hand column. As would 
be understood by a person skilled in the art, SOC is provided with 
CM/ECM data 12, derived from the HRR data. 

Iterations through CM 12 and ECM 14 are run with small loop 98 

15 replacing estimated 0,^, Q^, with Qf;, Q^^ converging to Q^, as 

taught above. An example of after conviction of a small iteration 

loop is 1661.9 kg/hr. An example of the resulting data file provided to IM 
10 for a subsequent big loop 162 is: 
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Symbol 


Value 


KrM 

ertg 


1 /M-u.u rpm 


'^brake 


1493.4 Nm 


^brake 


281.5 kW 


BSFC 


0.2177 kg/kWh 




3.05 bar 




441.8K 


Qin 


1661.9 kg/h 






SOI 


-30.0 CA 


Qegr 


415.0 kg/hr 


T„ 


818.0 K 


Pa 


2.855 bar 



Not all of data is used by IM 10 during each loop as described 
above. The above is also indicative of the data generated by ECM should 
the iteration of big loop end. 
5 ECM 14 with HRR and SOC provides i^^, an example of which 

is: 



-150 


3.02653 


-149.5 


3.03322 


-149 


3.04004 


-148.5 


3.047 


-148 


3.0541 






567 


2.99286 


567.5 


2.99874 


568 


3.00475 


568.5 


3.01089 


569 


3.01716 


569.5 


3.02357 



Where CA is provided in the first column and P^^ , in bar, is 

10 provided in the second column. Engine performance data, engine output 
data 38, ECM/EM data 33 and IM/ECM/EM data 25, is provided where no 
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further iterations of big loop are required as determined by €p after a 
comparison of P^f with P^,"* , an example of which is: 



Svmbol 


Value 


RPM 

^^^^ eng 


1740.0 rpm 


^ brake 


1493.4 Nm 


^ brake 


281.5 kW 


BSFC 


0.2177 kg/kW 


p 

^ CO 


3.05 bar 




441.8 K 


a. 


1661.9kg/h 


Q/uel 


51.500 kg/h 


SOI 


-30.0 CA 


Qegr 


415.0 kg/hr 


^' 


818.0 K 




2.855 bar 




2.969 bar 



5 After iterations through big loop 1 1 8, an example of P^, fi^om 

ECM 14, when P^, is set to P^, as described above, is: 



-360.000 


1 


.749 


-359.500 


1 


.751 


-359.000 


1 


.756 


-358.500 


1 


.762 


-358.000 


1 


.770 


-357.500 


1 


.778 


-357.000 


1.786 


357.000 


1. 


784 


357 .500 


1 . 


773 


358.000 


1. 


764 


358.500 


1 


757 


359.000 


X 


752 


359.500 


1 


749 


360.000 


1 


749 



Wifli the con^lete set of IM/ECM/EM data 25, ECM/EM data 16 
10 and EM system data 39, emissions data 41 is generated from EM 16. An 



wo 2005/028837 



PCT/CA2004/001735 



example of JVO, based on the above is 295.3 ppm. 

A higher number of iterations through big loop 118 helps to refine 
the accuracy prior to full generation of the ultimate output data, engine 
output data 38 and emissions data 41, from ECM 14 and EM 16. As well, 
5 the basis of this output data, data from IM 10 and CM 12, IM/ECM/EM 
data 25, IM/CM data 27 and CM/ECM data 29, is refined as well. 
Typically four iterations of big loop for an example such as the one above 
would be reasonable. 

CM 12 developed for directly injected gaseous fuel combustion is 
10 useful for the subject invention. The desired pressure at which gaseous 
fuels are injected into the combustion chamber of an internal combustion 
engine are dependent on P^^ to a far greater extent than is the case for 

liquid direct injection engines. Diesel fiiel is injected at pressures that can 
be an order of magnitude higher than is the case for natural gas. As such, 
15 the range in /J.^; is not as important to IM 10 accuracy. Therefore, when 

considering a gaseous-fuelled engine, lEM 8 incorporating P^^ derived 

from ECM 14 after completion of big loop 118, allows improved accuracy 
in IM 10 as well as improved functionality in IBM 8. Where liquid fuels 
are used however, there is less need for refining the value of P^^ as is 

2 0 done through the iteration described above. As such, there has been less 

need to deal with this aspect of liquid direct injection engines. 

For the purposes of this application, the sub-modules, IM 10, CM 
12, ECM 14 and EM 16 include the prior art examples discussed above 
and may vary depending on the application of the user: see background 
25 discussion. 

While particular elements, embodiments and applications of the 
present invention have been shown and described, it will be understood, 
that the invention is not limited thereto since modifications may be made 
by those skilled in the art without departing from the scope of the present 

3 0 disclosure, particularly in light of the foregoing teachings. 
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Nomenclature Table 

Definition 



Abbreviatio 
n 

ffiM 

IM 

CM 

ECM 

EM 

Q'o. 



GEO^ 

pn 
P 

SOI 
RPM 



and 



eng 



amb 

GEO^ 



eng 



amh 
cgr 

VGT^ 
p 

'■brake 
"^hmke 

BSFC 
NO, 

PM 

Qin 



Integrated engine model 
Injector module 
Combustion module 
Engine cycle module 
Emissions module 
Initial cylinder pressure 
Initial intake flow rate 
Initial EGR rate 
Injector geometry 

Iterative cylinder pressure 
Rail pressure 

Commanded or desired start of injection 

Commanded or desired pulse width 

Speed of engine 

Ambient pressure 
Engine geometry 

Ambient temperature 
Effective EGR valve flow area 

Variable geometry turbine position 

Brake power 
Brake torque 

Brake specific fuel consunq)tion 
Emissions of oxides of nitrogen 

Particulate matter emission 
Intake flow rate 
Cylinder pressure 
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Abbreviatio Definition 
n 

Turbine inlet temperature 

P^^ Turbine inlet pressure 

T^^ Compressor outlet temperature 

P^^ Compressor outlet pressure 

Turbine speed 

O EGR flow rate 

Q^^^ Exhaust flow rate 

HRR Heat release rate 

ROI Rate of injection 

SOI Start of injection 

Qfuei Fuel supply rate 

SOC Start of combustion 

Qf^ Iterative intake flow rate 

Iterative EGR flow rate 

/J^ Intake valve closing time cylinder pressure 

T^^ Intake valve closing time cylinder temperature 

Q^^^ Previous cycle intake flow rate 

Q^*' Previous cycle EGR flow rate 

Sq Intake air flow rate error 

€q^^ EGR flow rate error 

Sp Cylinder pressure error 

PW Pulse width 

P^^ Intake manifold pressure 

Effective volvune of cylinder 

Exhaust tenq)erature 
T^^ Temperature EGR 

Ar Flow restriction of intake and exhaust 
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Abbreviatio 
n 

X 

rpdmd 
nntax 

PM 



dmd 



NO 



cmd 

SOI 



,dmd 



cmd 



PW. 



cmd 



rail 

vgt;„ 



* egr 

CONP 
CONP'''^ 



SOI 



cmd 



PW. 



cmd 



^ rail 

vgt: 



pes 



A" 



SOL 



eng 



peng 
*rail 



Definition 

Air to fuel ratio 

Demanded conq>ressor outlet temperature limit 
Demanded turbine inlet temperature limit 
Maximum cylinder pressure 

Demanded particulate matter emission 

Demanded emissions of oxides of nitrogen 

Commanded brake torque 

Initial commanded start of injection 

Initial commanded pulse width of injection of 

fiiel 

Initial commanded rail pressure 

Initial commanded variable geometry turbine 

position 

Initial commanded effective EGR valve flow area 
Control parameter 
Demanded control parameter 
Iterative commanded start of injection 
Iterative commanded pulse width of injection of 
fuel 

Iterative commanded rail pressure 

Iterative commanded variable geometry turbine 

position 

Iterative commanded effective EGR valve flow 
area 

Start of injection control command to engine 
Pulse width of injection of fuel command to 
engine 

Rail pressure command to engine 
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Abbreviatio Definition 
n 

VGT^f Variable geometry turbine position command to 

engine 

A^^ Effective EGR valve flow area command to 

engine 

SOI'^ Final iterative commanded start of injection 

PJVj^ Final iterative commanded pulse width of 

injection of fuel 
p^^l Final iterative commanded rail pressure 

VGTp^l Final iterative commanded variable geometry 

turbine position 



A^'"' Final iterative commanded effective EGR valve 

flow area 



